The extremely red host galaxy of GRB 080207 
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ABSTRACT 

We present optical, near-infrared, and Spitzer IRAC and MIPS observations of the host galaxy of 
the dark gamma-ray burst GRB 080207. The host is faint, with extremely red optical-infrared colors 
(R - K = 6.3, 24//m/7?-band flux ~ 1000) making it an extremely red object (ERO) and a dust-obscured 
galaxy (DOG). The spectral energy distribution (SED) shows the clear signature of the 1 .6 p.m photometric 
"bump", typical of evolved stellar populations. We use this bump to establish the photometric redshift Zphot 
as 2.2^Q 3, using a vast library of SED templates, including M 82. The star-formation rate (SFR) inferred 
from the SED fitting is ~119M0yr"\ the stellar mass 3 x 10'' Mq, and Ay extinction from 1-2 mag. 
The ERO and DOG nature of the host galaxy of the dark GRB 080207 may be emblematic of a distinct 
class of dark GRB hosts, with high SFRs, evolved and metal-rich stellar populations, and significant dust 
extinction within the host galaxy. ^ 

Subject headings: galaxies: high-redshift — galaxies: ISM — (ISM:) dust, extinction 



1. Introduction 

Gamma-Ray Bursts (GRBs) without an optical 
counterpart of their X-ray afterglow are usually called 
"dark bursts"; a more precise definition based on their 
X-ray to optical spectral energy distribution (SEP) has 
been proposed ( Jakobss on et al ] |2004IRo1 et"al]|2005[ 



Ivan der Horst et al. 2009h . It is clear now that in most 
cases extinction by dust in their host galaxies makes 
these afterglows optically dim tareiner et al. 201 1), 
while cosmological Lyman drop out (high redshift) 
or in trinsic faintness are the exception rather than the 
rule dCenko et al.ll2009l: IPerley et al.lEo09llRossi etaT 



2011 ) 



The long-duration GRB 080207 is one of these 
truly dark GRBs. It was detected by Swift/BAI on 
2008 Februai-y 7 at 21:30:21 UT. SwiftPiRT began ob- 
serving the field 124 s after the BAT trigger and found 
a bright X-ray afterglow, with a po sitional uncertainty 
radius of I'.'A (iRacusin et al.ll2008a) . No optical/near- 
infrared afterglow was detected for GRB 080207, de- 
spite heroic efforts by several ground-based observa- 
tories. The deepest limit relative to the fading X-ray 
afterglow was achieved by the Zeiss-600 telescope at 
Mt.Terskol observatory which did not detect the after- 
glow down to j^AB > 20.5 at 1.69 hr after th e burst 
jAndreev et ai1l2008l) . iRacusin et all (l2008blcl) report 
that the XRT light curve between 1.30hr and 4.72 hr 
after the GRB start time declines monotonically with 
a power-law of index ax - 1.85 + 0.10. Therefore, we 
have assumed that the optical light curve in the same 
time interval is also decreasing monotonically, so that 
its behavior is s u fficien tly regular that the criterion of 
Jakobsson et al.l (2004) to define GRB "darkness" - 



based on the comparison of the X-ray and optical flux 
levels at 11 hr - can be applied to the optical upper 
limit Rab > 20.5 measured at 1.69 hr after the trigger 
The X-ray flux at 1.69 hr is -0.35 c/s (~ 3.1 x 10"" 
ergs ' cm^), which, assuming the spectral index fitted 
by Racusin et al. to the average XRT spectrum in the 
above time interval, /3x - 1 .4 + 0.2, corresponds to 
a 1 keV flux of 4.2 fi}y. Together with the simulta- 
neous optical upper limit, this yields an optical-to-X- 
ray index of fipx < 0.27, which leads to a dark GRB 
classification, according tolJakobsson et al.l (120041) and 
van der Horst et al.1 (l2009l) . 



In the vicinity of the XRT error circle, iRossi et al 



(1201 lb found two very faint host candidates, one bet- 
ter visible in VIMOS//?-band and the second brighter 
in K. However, the Chandra satellite observed the 



field of GRB ( )80207 nine days afte r the burst detec- 
tion by Swift jRacusin et al. 2008a ). and was able to 
localize the GR B X-ray afterglow with very high ac- 
curacy of ^.'67 (lEvans etal.ll2010l) : a = 13:50:02.97, 
5 = 07:30:07.8 (J2000). This posi tion coincid e s with 
one of the candidates found by iRossi et alj (1201 ll) 
(source "B"): a very faint (Rab ~ 26.5 mag), very red 
[R-K >6, (R - K)ab ~ 4.7] galaxy. The other candi 



date is bluer (dubbed "A" by Rossi et aD f201 1). with 
(R - K)ab ~ 2.1, and located just outside the 90% 
XRT error circle, ~2'.'3 northeast of the Chandra po- 
sition. Hence we associate the first source (B) with 
GRB 080207. 

In this Letter, we present proprietary and archival 
optical and near-infrared (NIR) observations, com- 
bined with archival Spitzer IRAC and MIPS data of 
the host of GRB 080207. We estimate the photometric 
redshift Zphot of the host by fitting the obs erved SEP 



with a vast library of GRASIL r nodels (|Silva et al 
'l998'; Iglesias-Paramo et al. 2007 ; Michalowski et al 



i2008, 2010). Pespite its optical faintness, because 
of the clear photospheric NIR bump at rest-frame 
/i=1.6/im, we are able to estimate the redshift of the 
host of GRB 080207 with a precision of +0.3. Sec- 
tion |2] describes the data reduction and the methods 
used to derive the photometry. The SEP models are 
presented in § [3] and we discuss the results and their 
implications in § HI Throughout the paper, we assume 
a Q,„ = 0.3, Qa - 0.7 cosmology, with Hubble constant 
//() = 70kms"' Mpc"'. 

2. Data and photometry 

In the course of an analysis of GRB host galaxy 
SEPs (Hunt et al., in preparation), we culled the 
Spitzer archive for observations of GRB 080207. It 
was observed by A. Levan (PIP 50562) more than 
one year after the burst in all four IRAC channels 



(IFazio et alj 120041) an d about five mon ths after the 
burst at MIPS 24;um (Rieke et al."2004). R- and K- 
band data were taken from .Rossi et al.. (.201 1.) . B- and 
7? -band images were acquired at the Large Binocular 
Telescope (LBT) in the course of our approved ob- 
serving program; /' and z' images were taken from 
the Gemini archive, and an i/-band image from the 
HST/NICMOS archive. Finally, we retrieved the 
ESO/SINFONI NIR integral field spectral data from 
the ESO archive, which provided a broad-band J con- 
tinuum point. Altogether we analyzed 12 photometric 
data points, from 0.4 jum to 24 fim, in order to establish 
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the photometric redshift and the general properties of 
the host of GRB 080207. The magnitudes given be- 
low are corrected for foreground Galactic extinction 
assuming E(B - V) = 0.023 mag dSchlegel et al.ll 19981) 
and a ratio of total-to-selective extinction of Ry - 3.1. 

2.1. Optical and near-infrared observations 

We observed the field 2 years after the burst with 
LBT/LBC in the B and R bands, and with VLT/VIMOS 
in R and VLT/ISAAC in K, bands. Optical and NIR 
data were reduced using standard IRAF taskfl and 
analyzed through aperture photometry. We used an 
aperture size twice the Full Width Half Maximum 
(FWHM) of the stellar PSF Optical VLT and LBT 
images were calibrated using Landolt stars in the 
field PG1047-H3. ISAAC fields were calibrated using 
2MASS field stars. 

Within the Swift/XRT error circle of 1'.'4, and coin- 
cident with afterglow localized by the Chandra satel- 
lite, an extended object (r.'6x0'.'9) is visible in our 
LBC, VIMOS 7;-band, and ISAAC K,-hnnd images; 
it is thus identifie d as the optical co unterpart of the 
GRB host galaxy ( IRossi et al.ll201 lb . The galaxy is 
marginally detected in both LBC and VIMOS R band 
with Rab = 26.43 + 0.37, but is clearly visible in our 
ISAAC images (see Fig. [1] Ks,ab = 21.77+0.14). 
The object is absent in our B-band LBT image, which 
is however compromised by worse seeing. 

The host galaxy is also clearly visible in an archival 
Gemini/GMOS z' image obtained 6 months after the 
burst during very good sky conditions, but not in the 
i' image taken immediately after the GRB. We obtain 
i'j^j^ > 24.6 and z^^g = 25.02 + 0.02. 

We reduced archival J (30 min exposure) and H + K 
(45 min) band SINFONI integral field spectra around 
the host galaxy, following all the typical reduction 
steps applied to near-IR spectra. After background 
subtraction, the data were flat-fielded, corrected for 
dead/hot pixels, wavelength and flux calibrated using 
telluric stars. The resulting datacubes were collapsed 
to obtain broadband continuum images; although the 
final S/N in the H + K band was too low, the derived 
J band image was analyzed using standard aperture 
photometry to extract the flux. The spectra themselves 
were of insuflicient sensitivity for reliable detection of 



[On], [Om], or Ha emission lines (see ®. 

We reduced the HST/NICMOS F160W (NIC3) im- 
age starting with the standard flat-fielded calibrated 
frames after correcting for the pedestal effect. These 
frames were then "drizzled" onto the final mosaic 



dFruchter & Hook 2002 ). using a pixel filling factor of 
0.65, and maintaining the original pixel size of 0'.'2. 
The host galaxy of GRB 080207 is clearly visible in 
this image also, and using the same methods as de- 
scribed above, we obtain a flux of 2.9+0.7 juJy (see Ta- 
ble[B. 

2.2. Spitzer IRAC and MIPS observations 

For both IRAC and MIPS observations, we started 
the data reduction with the Basic Calibrated Data 
(BCD) and the corresponding masks (DCE masks), 
furnished by the Spitzer Science Center (SSC) pipeline. 
The image mosaicing and source extraction package 
(MOPEX, Makovoz & Marleau 2005) was used to 
co-add BCD frames for each source, and mask out 
bad pixels. For MIPS-24, we discarded the first two 
frames of each observation sequence, as advocated 
by the MIPS Instrument Handbook. Also for MIPS, 
we incorporated the sigma-weighting algorithm since 
it gave less noisy MIPS mosaics than without. The 
frames were corrected for geometrical distortion and 
projected onto a fiducial coordinate system with pixel 
sizes of r.'20, roughly equivalent to the original IRAC 
pixels. The quality of our IRAC and MIPS maps is 
substantially better than that provided by the Spitzer 
archive. 

Fi gure [H show s the IS AAC TiT-band image taken 
from 



Rossi et al 



'iRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 



(120111. together with the IRAC 
3.6 jum and MIPS 2Apm images presented here. The 
Chandra afterglow position is marked with an x. The 
brightest IRAC pixels correspond to the Chandra po - 
sition, and to source "B" found bv lRossi et alj (1201 lb . 

In all IRAC and MIPS Spitzer images we detected 
only one object within the Chandra error box which 
we identify as the GRB 080207 host galaxy. The ob- 
ject is not spatially resolved in the Spitzer images. We 
estimated the flux density in all filters using a small 
circular aperture with a radius of 2 pixels in IRAC 
and 3 pixels in MIPS; these measurements were per- 
formed using the phot task in IRAF. Aperture correc- 
tions have been applied to account for the extended 
size of the PSF relative to the small apertures. 

Due to the possible contamination from nearby ob- 
jects (see Fig. [U, we checked the resulting fluxes 
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Fig. 1.— From left to right: the host galaxy of GRB 080207 in K band (^ Rossi et alJboilb . in IRAC/3.6//m, and 
MIPS/24 yum. North is up. East is left; the images are 30 " on a side. The XRT error circle is shown, and the Chandra 
position is marked with x. 



by using PSF fitting with the DAOPHOT software 



ins 
2d . 



(lStetsonlll987.) to perform the source extraction and 



photometry. The results from the two methods are per- 
fectly consistent. Table [T] reports the fluxes we have 
derived and use in the SED fitting described below. 

3. Modeling the SED 

Wifli {R - K)ab - 4.7 [{R - /:)vega ~ 6.3], the host 
galaxy is very red, among the reddest host galaxies of 
a LGRB ever observed. Such red R - K colors clas- 
sify the host galax y of GRB 080207 a s an extremely- 



red object (ERO: lElston et alj|1988l) . The R-24nm 
color is also sufficiently red (24 fim/R flux ratio ~994) 
to classify it as a dust-obscured ga laxy (DOG: 24 fim/R 
flux ratio > lOOO. IPey e"taDl2008h . 



Given the extremely dim optical fluxes, the strongest 
constraint on the SED we present here is the spectral 
"bump" at rest-frame 1.6yum (see Fig. |2|l. This bump 
is due to the minimum in the opacity in the ion 
found in the photospheres of cool stars, and found in 
virtually all stellar populations older than ~100Myr; 
it provides a conspicuous signa ture with which to con- 
strain photometric r edshifts (ISimpson & Eisenhardt 
19991: ISawickil l2002h . The host of GRB 080207 is 
one of the few GRB hosts to show this feature (an- 
oflier one is ESQ 184-G 82, the host of GRB 980425, 
Michalowski et al. |2009 >). Because of its extremely 



red colors, the host of GRB 080207, is either a pas- 
sive early-type galaxy at redshift 1 < z < 2, or a 
dusty starburst at higher redshift. Because LGRBs 
occur in star-forming gala xies (e.g., Le Floc'h et al] 
20031: ISavaglio et alJ l2009h . we concentrated on the 



second alternative, and considered a library of - 
35 000 GRASIL star-forming galaxy templates, in 
eluding M 82 and Arp220 and oflier GRB host tem- 



2007 



plates fSilva et al. 1998; Iglesia s-Paramo et al 
iMichalowski et al. , 2008. 2010.) . 

SED fitting in vLy space was performed by renor- 
malizing the templates, to minimize over the ob- 
ser ved SEeB This was done by redshifting the models 
over a range of redshifts from to 6, in increments of 
0.1. Lyman fo rest attenuation was treated according to 



Mada ul (fl995 ). We used t he fitti n g method described 



Iglesias- Paramo et alj ( l2007h : iMichalowski et al. 



20081 l20ld) and one developed independently by 
our group; the two approaches give virtually identi- 
cal results. The best fit was obtained for a starburst 
template, with roughly the same parameters as M 82, 
at a redshift z - 2.2 (mean residuals 0.1 yuJy, reduced 
xl = 0.3). The best-fitting SED is shown in Fig. |2] 
and the trends with redshift of xl for representative 
GRASIL templates are shown in Fig. [3] After a statis- 
tical analysis of the x^ distributions, we conclude that 
the best estimate of the photometric redshift for the 
host of GRB 080207 is z = 2.2+jJ-2 (errors are Icr, and 
were calculated including the choice of template as a 
free parameter). 

The most salient feature of our SED fit is its ability 
to well fit the MIPS point, the B and /' upper limits, and 
the R, z', and 7-band fluxes. Even though the 24fj.m 
data point is not included in our fitting procedure, it is 



^In order to better weight the NIR bump in the fits, we excluded both 
the 8/jm IRAC point and the 24 /jm MIPS point in the fits; the fits 
were more stable with this exclusion. 
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Table 1 

Photometry of the host galaxy of GRB 080207 





w 


Flux" 


Uncertainty 


Telescope/ 


Filter 




(/jJy) 




Instrument 


B 


0.44 


<0.I5'' 




lbt/lbc 


R 


0.66 


0.093 


0.026 


vlt/vimos 


/' 


0.78 


<0.32'' 




Gemini/GMOS 




0.925 


0.35 


0.06 


Gemini/GMOS 


J 


1.2 


1.6 


0.3 


VLT/SINFONI 


H 


1.60 


2.9 


0.7 


hst/nicmos 




2.2 


7.3 


1.0 


vlt/isaac 




3.550 


14.40 


0.31 


Spitzer/IRAC 




4.493 


15.51 


0.44 


Spitzer/IRAC 




5.731 


18.53 


1.58 


Spitzer/IRAC 




7.872 


12.52 


1.76 


Spitzer/IRAC 




23.680 


92.43 


6.50 


Spitzer/MIPS 



"Blueward of 3.5 ^m, these values have been conected for Galac- 
tic extinction as described in the text. 

^30" upper limits. 




Observed X [fj.m] Photometric redshift 



Fig. 2. — The observed SED (corrected for Galactic 
extinction as described in the text) of the host galaxy 
of GRB 080207 plotted against observed wavelength. 



Fig. 3.— GRB 080207 host galaxy xl of the NIR 
bump plotted against redshift. The best-fit redshift 
z = 2.2 is indicated by a dotted line. 
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well predicted by the best-fit SED. At z = 2.2, MIPS- 
24 falls roughly on the redshifted polycyclic aromatic 
hydrocarbon (PAH) feature, ubiquitous in metal-rich 
starbursts. The faint optical fluxes are clearly subject 
to significant extinction; the GRASIL starburst mod- 
els predict Av~l mag. However, this could be a lower 
limit since these opacities are integrated over the en- 
tire galaxy in the models, and consider the opacity of 
the molecular clouds, together with their filling factor 
over the galaxy disk. 

4. Discussion and conclusions 

We can use the GRASIL templates to infer the total 
infrared (IR) luminosity (Lir) the star-formation rate 
(SFR), and the steUar mass (M*). From the SED fitting 
procedure, we obtain a total Lir~ 6.8 x 10" + 25% Lq, 
and SFR ~ 119+25% Moyr^^ This is roughly an 
order of magnitude lower than the SFRs implied by 
the few detections of GRB hosts in the submillime- 
ter( Tanvir et al but would definitely quahfy the 

host of GRB 080207 as a dusty luminous star-forming 
galaxy. 

In the absence of extinction, the inferred SFR is suf- 
ficiently high that both [O ii] and Ha would have been 
detected at > 5cr with the SINFONI spectra. The (con- 
servative) 5cr upper limits we obtain (in J for [On] 
and K for Ha) give a SFR- 30 - 40Mo yr ^' (using the 



SFR calibrations of ISavaglio et al.l [20091) . taking into 
account the uncertainty in redshift, and the impact of 
the sky lines on the noise in the observed spectrum. 
This is a factor of 3 to 4 times lower than the unattenu- 
ated SFR inferred from Lir of the SED analysis; if the 
latter value is the true one, the extinction of the emis- 
sion lines would be A v> 1-2 mag, consistent with the 
SED fitting. 

There is some indication that the extinction could 
be even higher around the GRB 080207 location. As- 
suming there is a spectral break between the opti- 
cal and X-ray band due to electron cooling (vc), the 
optical spectral slope soon after the trigger will be 
/3opi - Px - 0.5, and the R magnitude is fainter for 
higher values of the frequency break. If vc coincides 
with the lower boundary of the XRT spectral range (0.3 
keV), the extrapolated R band magnitude is ~ 15.4, 5. 1 
magnitudes brighter than the photometric upper limit 
(7?AB > 20.5). The optical flux deficit is thus signifi- 
cant and, if due to dust extinction, it supports indepen- 
dently the results we have obtained from the analysis 
of the galaxy SED. 



The implied luminosity and SFR places the host 
galaxy of GRB 080207 into a subgroup of GRB host 
galaxies characterized by violent star formation and 
clearly distinct from the majority of the GRB host 



galaxy population (c.f., Svensson et al. 2010l) . A 



representative example is the host of GRB 010222, 
which is supposed to have Lboi ~ 4 x lO^^L p 
and a SFR of about 600 Moyr ' jFrail et al.1 [20021) . 
Also the submm-detected hosts of GRBs 980329, 
980703, 000210 and 000418 belong to this c ategory 



(IBerger et all 120031: iMichalowski et all l2008h . while 



the host of the first dark burst, GRB 970828, only 
has a comparatively modest SFR of 24;^|^M0yr"' 
dLe Floc'h et alJl2006h . 

The steUar mass, M., of the host of GRB 080207 
inferred from the SED fitting is ~ 3.2 x 10" +25% Mg, 
one to two orders of magnitude higher tha n those 
found for GRB hosts at z <1.5 (e.g.. Savaglio et al. 
2009h . The SED fit gives an absolute rest-frame 



K magnitude of the GRB 080207 host: M/r(AB) ~ 
-25.2. Fr om the correla ti on be tween M, and Mk 
derived by ISavagho eta D (l2009l) . we infer M, ~ 
4.0 X 1O"M0, consistent within the uncertainties of 
the value derived from the SED analysis. 

Figure |4] shows stellar masses plotted against red- 
shift for GRB host gala xies and a set of GOODS galax- 
ies (taken from Savagl io et alj|2009.) . As can be seen 
in the figure, the mass of the host of GRB 080207 is 
comparable to the most massive GRB host galaxies 



ever measured (GRBs 02012 7, 080325: iBerger et al, 



20071; [Hashimoto et all I2OI0I) . Besides being quite 



massive, the hosts of bofli GRB 020127 and 080325 



are al so EROSs (IBerger et al.. .2007.; Hashimoto et al, 



201 oh . similar to the host of GRB 080207. The SFR 



and stellar mass of the GRB 080207 host are consis- 
tent with the correlation found for star-fo rming BzK- 
selecte d galaxies at z = 2 discussed by iDaddi et al. 



(I2007h . The specific SFR, SSFR, of the GRB 080207 
host is comparable to other galaxy populations aiz — 
2, with SSFR ~ 0.4 Gyr '. 

The metal abundance i mplied from the mass- 
metallicity relation at z = 2 dErb et al.l2006l ; Havashi et al 
2009) places the massive host of GRB 080207 at 
roughly solar metal abundance. This galaxy is more 
chemically evolved than typical GRB hosts which (at 
least at lower reds hifts) tend to be me tal poor and 
much less massive ( Savaglio et al. 20091) . Moreover, 
the GRB 080207 galaxy seems to be older than many 
GRB hosts. The mass of the stars in the burst is <3% 
of the total stellar mass in the galaxy. This is consis- 



6 



12 - 



1 



O 




■ GRB EROs 

•Other GRB hosts 

'• Other dark GRB hosts 

• GOODS South galaxies 

• GRB 080207 host 



redshift 

Fig. 4. — Stellar mass vs. redshift for GRB hosts 
and other galaxy populations. Data are taken from 
fSavaglio et al..(.2009.) . 



tent with the clear presence of the NIR bump in the 
SED, and implies that the mean stellar age is relatively 
old, ~l-2Gyr, unlike man y GRB hosts which tend 
to be young, ^lOMy r (e.g.. iMichalowski et al.ll2008 : 



Levesque et 



, < lUMy r 
al] l2010l) . 



With its extremely red R-K color, (R-K)yega ~ 6.3, 
the host of GRB 080207 is the most extreme ERO 
found to date for a GRB host galaxjH Moreover, 
the unusually high ratio of 24//m-to-7?-band flux of 
the GRB 080207 host makes it also a dust-obscured 



galaxy, or DOG dDev et al.1 120081) . To our knowl- 
edge, this is the first GRB host galaxy to be so classi- 
fied. Such red host colors at z >2 stand out from the 
typical sub-luminous blue galaxies tha t tend to dom 



inate lower-z GRB host populations dSavaglio et al 



20091) . In fact the colors and other properties of 



the host of GRB 080207 are more typical of dusty 
z -2 galaxies de tected in the IR and submillimeter 



(IPope et al 



ies (e.g., 
or Distant 



2008h. star - forming feA"-selected galax - 



Daddi et al.1 120051: iGrazian et al. 2007), 



Red Galaxies (e.g., Ivan Dokkum et al 



20041) . Indeed, there is more and more evidence (e.g. 



GRBs 20127, 080325 080 607: Berger et al.M2007l: 
[Hashimoto et all I2OI0I: IChen et al.. .2010.) that some 
dark bursts, those with extremely red optical-IR col- 
ors, may point to a distinct GRB host population: one 
that harbors more evolved, more metal enriched stellar 



^The host of GRB 020127 has (R - ^:)vega ~ 6.2 (at z = 
iBerger et alj|2007l). and the host of GRB 0301 15 (R - S:)vega 
(z = 2.6: lLevan et alj200^ . 
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populations, but with high SFRs and significant dust 
extinction. 
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